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The issues of quintessence and cosmic acceleration can be discussed in the 
framework of higher order curvature and torsion theories of gravity. We can 
define effective pressure and energy density directly connected to the curvature 
or to the torsion fields and then ask for the conditions to get an accelerated 
expansion. Exact accelerated expanding solutions can be achieved for several 
fourth order curvature or torsion theories so that we obtain an alternative 
scheme to the standard quintessence scalar field, minimally coupled to grav- 
ity, usually adopted. We discuss also conformal transformations in order to 
see the links of quintessence between the Jordan and Einstein frames. Fur- 
thermore, we take into account a torsion fluid whose effects become relevant 
at large scale. Specifically, we investigate a model where a totally antisym- 
metric torsion field is taken into account discussing the conditions to obtain 
quintessence. We obtain exact solutions also in this case where dust dominated 
Friedmann behavior is recovered as soon as torsion effects are not relevant. 
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1 Introduction 



One of the astonishing recent result in cosmology is the fact that the universe is acceler- 
ating instead of decelerating along the scheme of standard Friedmann model as everyone 
has learned in textbooks. Type la supernovae (SNe la) allow to determine cosmologi- 
cal parameters probing the today values of the Hubble constant Hq and the deceleration 
parameter go El- Besides, data coming from clusters of galaxies at low red shift (includ- 
ing the mass - to - light methods, baryon fraction and abundance evolution) 0, and 
data coming from the CMBR investigation (e.g. BOOMERANG)]^ give observational 
constraints from which we deduce the picture of a spatially flat, low density universe 
dominated by some kind of non-clustered dark energy. Such an energy, which is supposed 
to have dynamics, should be the origin of the cosmic acceleration. 
In terms of density parameter, we have 

^n^a^^er) ^ 0.3 , ^ 0.7 , fifc ~ 0.0 (1) 

where the matter is the non-relativistic baryonic and non-baryonic (dark) matter, A is 
the dark energy (cosmological constant, quintessence,..), k is the curvature parameter of 
Friedmann-Robertson- Walker (FRW) metric of the form 



ds^ = dt^ - a{tf 



1 — kr 



(2) 



where a{t) is the scale factor of the universe. 

The luminosity distance can be deduced from SNe la used as standard candles. 

For 2 < 1, by the luminosity distance dL — Hq^[z + (1 — go)-2^/2], the observational 
results indicate 

- 1 < go < (3) 

which is a clear indication for the acceleration. 

The deceleration parameter can be given in terms of density parameter and then we 
have for FRW models 

go = -T^ = ^(37 + l)^^(ma«er) - ^^A , (4) 

where 7 is the constant of state equation p = 7p. Immediately, we realize that acceleration 
or deceleration also depends on the value of 7. For standard fluid matter, it is defined the 
Zeldovich range < ^(matter) < 1 where ^{matter) = indicates dust (i.e. non-relativistic 
matter), ^{matter) = 1/3 radiation (i.e. relativistic matter). Another interesting case 
has been widely considered in literature, it is 7 = —1 which points out a scalar field 
fluid where dynamics is dominated by self-interaction potential or cosmological constant. 
Inserting a standard matter fluid into Friedmann-Einstein cosmological equations gives 
rise to decelerated dynamics. Due to this fact, non-standard forms of matter energy have 
to be taken into account if one wants to insert observations in a theoretical frame. 
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Several approaches can be pursued in order to realize this goal. All of them can be 
summarized into three great families: the cosmological constant, the variable cosmo log- 
ical constant and quintessence. Essentially they are linked but, in order to match the 
observations, several issue have to be satisfied. Below we give a short summary of this 
three pictures. 

The cosmological constant has become one of the main issue of modern physics since 
by fixing its value should provide the gravity vacuum state |Q, should make to understand 
the mechanism which led the early universe to the today observed large scale structures 

0. 0, and to predict what will be the fate of the whole universe (no-hair conjecture) 0. 
From the cosmological point of view, the main feature of inflationary models is the 

presence of a finite period during which the expansion is de Sitter (or quasi-de Sitter or 
power law): this fact implies that the expansion of the scale factor a{t) is superluminal (at 
least a{t) ~ t, in general a{t) ~ exp HqI where Hq is the Hubble parameter nearly constant 
for a finite period) with respect to the comoving proper time t. Such a situation arises 
in presence of an effective energy-momentum tensor which is approximately proportional 
(for a certain time) to the metric tensor and takes place in various gravitational theories: 

1. e. the Einstein gravity minimally coupled with a scalar field fourth or higher-order 
gravity g |, llOl In] scalar-tensor gravity pi pi. 



Several inflationary models are affected by the shortcoming of "fine tuning" that 
is inflationary phase proceeds from very special initial conditions, while a natural issue 
would be to get inflationary solutions as attractors for a large set of initial conditions. 
Furthermore, the same situation should be achieved also in the future: if a remnant of 
cosmological constant is today observed, the universe should evolve toward a final de 
Sitter stage. A more precise formulation of such a conjecture is possible for a restricted 



class of cosmological models, as discussed in |T^. We have to note that the conjecture 
holds when any ordinary matter field, satisfies the dominant and strong energy conditions 
T6| . However it is possible to find models which explicitly violate such conditions but 



satisfies no-hair theorem requests. Precisely, this fact happens if extended gravity theories 
are involved and matter is in the form of scalar fields, besides the ordinary perfect fluid 
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In any case, we need a time variation of cosmological constant to get successful in- 
flationary models, to be in agreement with observations, and to obtain a de Sitter stage 
toward the future. In other words, this means that cosmological constant has to acquire 
a great value in early epoch (de Sitter stage), has to undergo a phase transition with a 
graceful exit and has to result in a small remnant toward the future [0. The today ob- 
served accelerated cosmological behaviour should be the result of this dynamical process 
where the value of cosmological constant is not fixed exactly at zero. 

In this context, a fundamental issue is to select the classes of gravitational theories 
and the conditions which "naturally" allow to recover an effective time-dependent cos- 
mological constant without considering special initial data. 

The third approach is quintessence 0]. Quintessence is a time- varying, spatially 
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inhomogeneous component of cosmic density with negative pressure — 1 < 7q < 0. For- 
mally, vacuum energy density is quintessence in the limit 7q — — 1 so that the three 
approaches present in literature (cosmological constant, variable cosmological constant 
and quintessence) are strictly linked. 

However all of them claim for an ingredient which, a part a pure cosmological constant, 
comes from a matter-energy counter part. In this paper, we want to investigate if the 
quintessential scheme can be achieved in a geometrical way. 

There is no a priori reason to restrict the gravitational Lagrangian to a linear function 
of the Ricci scalar R minimally coupled with matter [2^. Additionally, we have to note 
that, recently, some authors have taken into serious consideration the idea that there are 
no "exact" laws of physics but that the Lagrangians of physical interactions are "stochas- 
tic" functions with the property that local gauge invariances (i.e. conservation laws) are 
well approximated in the low energy limit and physical constants can vary . This scheme 
was adopted in order to treat the quantization on curved spacetimes and the result was 
that the interactions among quantum scalar fields and background geometry or the grav- 



itational self-interactions yield corrective terms in the Einstein-Hilbert Lagrangian |21 



Futhermore, it has been realized that such corrective terms are inescapable if we want 
to obtain the effective action of quantum gravity on scales closed to the Planck length 
P^ . They are higher-order terms in curvature invariants as R^, R^^R^u, R^'^°''^R^iuai3, or 
nonminimally coupled terms between scalar fields and geometry as (p'^R. Terms of these 
kinds arise also in the effective Lagrangian of strings and Kaluza-Klein theories when the 
mechanism of dimensional reduction is working p3| . 

Beside fundamental physics motivations, all these theories have acquired a huge in- 
terest in cosmology due to the fact that they "naturally" exhibit inflationary behaviours 
and that the related cosmological models seem very realistic P, |12|. Furthermore, it is 
possible to show that, via conformal transformations, the higher-order and nonminimally 
coupled terms {Jordan frame) always corresponds to the Einstein gravity plus one or more 
than one minimally coupled scalar fields [Einstein frame) |^, ^ ^ so that these 
geometric contributions can always have a "matter" interpretation. 

As we will show below, quintessence can be achieved also in the framework of higher- 
order theories of gravity, that is in a geometrical way. 

Alternatively, taking into account torsion is a straightforward generalization to imple- 
ment concepts as spin in General Relativity I^S], However, it was soon evident that 
torsion, as considered e.g in Einstein-Cartan-Sciama-Kibble (ECKS) theory, does not give 
relevant effects in the observed astrophysical structures. Nevertheless it was found that for 
densities of the order of 1^)^'^ g/cm? for electrons and IQ^'^g/cm? for protons and neutrons, 
torsion could give observable consequences if all the spins of the particles are aligned. 
These huge densities can be reached only in the early universe so that cosmology is the 
only viable approach to test torsion effects However no relevant tests confirming 

the presence of torsion have been found until now and it is still an open debate if the 
space-time is Riemannian or not. Considering the cosmological point of view and, in 
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particular the primordial phase transitions and inflation ||38|, |39|, it seems very likely 



that, in some regions of the early universe, the presence of local magnetic flelds could have 
aligned the spins of particles. At very high densities, this effect could influence the evo- 
lution of primordial perturbations remaining as an imprint in today observed large scale 
structures. In other words, a main goal could be to select perturbation scales connected 
to the presence of torsion in early epochs which give today-observable cosmological effects 
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From another point of view, the presence of torsion could give observable effects with- 
out taking into account clustered matter. 

We want to investigate if the quintessential scheme can be achieved also by taking into 
account theories of gravity with torsion. 

2 Curvature Quintessence 

A generic fourth-order theory in four dimensions can be described by the action 

A= (fxyf^ \f{R) + ^(matter) 1 , (5) 



where f{R) is a function of Ricci scalar R and Li^natter) is the standard matter Lagrangian 
density. We are using physical units StiGn = c = h = 1. The fleld equations are 

f\R)Rap - \f{R)g.p = fiRy'^^g^.gp. - g^pg,.) + Ti7"^^^ , (6) 
which can be recast in the more expressive form 



n — T? ^ ^ r> — rp(curv) rp(matter) /«x 

— — -gafi-ti — -I- al3 + al3 i I ' J 



where 



and 



Ti^r^ = [f{R) - Rf'iR)] + f'm''''{ga,gp. - g.pg,.) \ (8) 

rp{matter) 1 rj^(matter) /q\ 

is the stress-energy tensor of matter where we have taken into account the nontrivial 
coupling to geometry. The prime means the derivative with respect to R. 

However, if f{R) = R + 2A, the standard second-order gravity is recovered. Reducing 
the action to a point-like, FRW one, we have to write 

A^curv) = / dtC{a, d; R, R) , (10) 
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(a a?' k 
R + 6{- + — + — 



where dot means derivative with respect to the cosmic time. The scale factor a and the 
Ricci scalar R are the canonical variables. This position could seem arbitrary since R 
depends on a, a, a, but it is generally used in canonical quantization |10|, The 
definition of R in terms of a, a, a introduces a constraint which eliminates second and 
higher order derivatives in action (|10D, and gives a system of second order differential 
equations in {a, R}. Action (|iy) can be written as 

A(curv) = 2vr2 j dt |aV(^) - A 

where the Lagrange multiplier A is derived by varying with respect to R. It is 

A = a'f'iR) . 

The point-like Lagrangian is then 

C = Ci^curv) + Ci^matter) = a' [f (R) - + 6ad^f'{R) + 

+ 6a^dRf"{R) - 6kaf'{R) + a^p(matter) , 



(11) 



(12) 



(13) 



where we have taken into account also the fluid matter contribution which is, essentially, 



a pressure term [|17 



The Euler-Lagrange equations are 



+ 



+ 



-P{tot) 



and 



fiR) 



„ fd ^ 

R + 6{- + — + — 



The dynamical system is completed by energy condition 

2 k 1 



(14) 



(15) 



Combining Eq. (|l4D and Eq.(|l6D, we obtain 

P{tot) + 3p(tot)J , (17) 
where it is clear that the accelerated or decelerated behaviour depends on the rhs. However 

P(tot) = P(curv) ~\~ P {matter) P{tot) = P{curv) + P{matter) , (18) 

where we have distinguished the curvature and matter contributions. 





1 r 


(!) 


~ ~6 - 
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From the curvature-stress-energy tensor, we can define a curvature pressure 

Picur.) = 77^ {2 (^) RfiR) + Rf"{R) + R'f"\R) - \ [f{R) - Rfm] , (19) 



f'{R) 

and a curvature density 

P{curv) ' 



1 rl 



- [f{R) - Rfm - 3 (-) RfiR)] . (20) 



From Eq.(|T^, the accelerated behaviour is achieved if 

P(tot) + ^P{tot) < , (21) 

which means 

1 

P{CUTV) > -^P(tot) , (22) 

assuming that all matter components have non-negative pressure. 

In other words, conditions to obtain acceleration depends on the relation 

Picur.) + Mcur.) = [r' f (R) + (^) RfiR) + Rf{R) ' ^ [f (R) - RfiR)]] , 

which has to be compared with matter contribution. However, it has to be 



(23) 



^^^^ = Kcurv) , -1 < Kcurv) < . (24) 

P{curv) 

The form of f{R) is the main ingredient to obtain this curvature quintessence. 

As simple choice in order to fit the above prescriptions, we ask for solutions of the 
form ^ 

f{R) = foR\ a(t)=ao(^) . (25) 

However, the interesting cases are for n 7^ 1 (Einstein gravity) and /? > 1 (accelerated 
behaviour). Inserting Eqs. ([25|) into the above dynamical system, we obtain the exact 
solutions 

p = 2; n = -l, 3/2; k = 0. (26) 
In both cases, the deceleration parameter is 

?o = , (27) 
in perfect agreement with the observational results. 
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The case n = 3/2 deserves further discussion. It is interesting in conformal transfor- 
mations from Jordan frame to Einstein frame ^ since it is possible to give exphcit 
form of scalar field potential. In fact, if 



ga(3 = f'{R)9al3 , 

we have the conformal equivalence of the Lagrangians 

R 1 



(28) 



-9 



foR'^' 



'-Vf,(fV^(f - Vo exp 




(29) 



in our physical units. This is the so-called Liouville field theory and it is one of the few 
cases where a fourth-order Lagrangian can be expressed, in the Einstein frame, in terms 
of elementary functions under a conformal transformation. It is possible to obtain the 
general cosmological solution which is 



a{t) = aQ[c4t^ + c^t^ + C2t^ + cit + Cq] 



1/2 



(30) 



The constants q are combinations of the initial conditions. Their values determine 
the type of cosmological evolution. For example, C4 7^ gives a power law inflation while, 
if the regime is dominated by the linear term in ci, we get a radiation-dominated stage. 



3 Torsion Quintessence 



In the ECSK theory, the affine connection is non-symmetric in its lower indices and the 
antisymmetric part 



i [be] 



Q a 



'be 



(31) 



is called torsion. Usually, we can divide such an antisymmetric part into three components: 
one of them is irreducible, while the other two can be set to zero |]42|, |43[|. This assumption 
yields the simplest theory containing torsion. 
Furthermore, we can express by a 4- vector 



a 



ubed Q 
£ J bed 



(32) 



the totally antisymmetric part of torsion. If one imposes to it the symmetries of a back- 
ground which is homogeneous and isotropic, it follows that, in comoving coordinates, only 
the component 0"° survives as a function depending only on cosmic time (see ^). 



For a perfect fluid, the Einstein-Friedmann cosmological equations can be written as 



a 1 / ~ 

- = -7 P + 3p 

a 



(33) 
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and 

'a\^ k 1 

(assuming SttG = c = 1) where a{t) is the cosmological scale factor and k = 0, ±1 is 
the spatial curvature constant. Energy density and pressure can be assumed in the forms 



-) +-=oP' (34) 



P = P + f\ P = P-f\ (35) 

where / is a function related to cr°, while p and p are the usual quantities of General 
Relativity. This choice can be pursued since we can define Sabc = S[abc] = f{t) where f{t) 
is a generic function of time which we consider as the source of torsion. For a detailed 



discussion of this point see [43 



As usual, we can define a stress-energy tensor of the form 

Tab^ = {P + P)UaUb - P9ab , (36) 

which, by Eqs. (^), can be splitted as 

rp(tot) rp(matter) rp{torsion) ('37\ 

^ab — -l-ab ~^ -'^ab ■ {"^ ' ) 

Due to the contracted Bianchi identity, we have 

tT^''' = ; (38) 
from which, we can assume that (cfr. ||45|| ) 

T,{matter);b rp{torsion);b ^ /'Qn^ 

ab -U, J-ab -U- [oy) 

In the FRW space-time, (|38|) becomes 

'p + 3H{p + p) = (40) 

which is 

p + 3H{p + p) = -2ff. (41) 
From Eqs.(|39D, both sides of (|4ll) vanish independently so that 

f{t) = fo = constant. (42) 

In other worlds, a torsion field gives rise to a constant energy density. Taking into account 
standard matter which equation of state is defined into the above Zeldovich range, we 
obtain 



8 



P = Po 



Oq 

a 



3(7+1) 



+ fo- 



(43) 



Inserting this result into the cosmological equations, we obtain, in any case, a mono- 
tonic expansion being /q > 0, > 0. The condition to obtain the accelerated behaviour 
is 

(44) 



p + 3p<2/2. 



so that acceleration depends on the torsion density. 
In a dust-dominated universe, we have 



P = Po 

and then the general solution is 

ait) = i 



ao 



+ fo 



P 



Jo ' 



(45) 



%Po\ 



1/3 



2/o 



cosh(/ot) - 1] 



(46) 



Obviously, if f^t 0, we have cosh(/ot) ^ (fo^Y and then a ~ t^/^, as it has to be. 



4 Conclusions 

In this paper, we have shown that the quintessence "paradigm" can be recovered in the 
framework of higher-order curvature and torsion theories of gravity. In other words, as it is 
possible for inflationary models, we can ask for a sort of curvature or torsion quintessence 
which can be recovered by taking into account curvature and torsion geometric invariants. 
The interest of this approach is that quintessence could be related to some effective theory 
of quantum gravity where such invariants are widely derived [HH . 



Besides, this results tell us that quintessence can be obtained without considering 
additional scalar fields in the dynamics but only assuming that the space-time is f/4 
instead of V4. A preminent role is played by the type of torsion we are going to consider. 
In our model, we assume that torsion density is a function of cosmic time only by imposing, 
as standard, the homogeneity and isotropy of background. However, such a density should 
be comparable to the observed limits of dark energy [i.e. Qa ~ 0.65 ^ 0.7) in order to 
give relevant effects. Another point is that such a torsion quintessence should match the 
issues of cosmic coincidence as scalar field quintessence. This point strictly depends 



on the type of torsion since torsion can be or not related to the spin density [Q. In 
the first case, the spin of baryonic and non-baryonic matter would rule the dark energy 
(torsion) density. In forthcoming investigations, we will face these problems. 

In any case, the models have to be improved by the comparison with observations in 
order to see if it is possible to constrain the form of f{R) and torsion field without ad 
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hoc choices. In this sense, some results are present in hterature where the form of f{R) 
is selected by the CMBR constraint p4 . 



ACKNOWLEDGEMENTS 

The authors are grateful to A. A. Starobinsky and S.A. Bludmann for the useful dis- 
cussions and comments on the topics. 

References 

[1] B.P. Schmidt et al. Ap. J. 507, 46 (1998). 
A.G. Riess et al. Ap. J. 116, 1009 (1998). 
S. Perlmutter et al. Ap. J. 483, 565 (1997). 
S. Perlmutter et al. Nature 391, 51 (1998). 
S. Perlmutter et al. Ap. J. 517, 565 (1999). 

[2] B. Chaboyer et al. Ap. J. 494, 96 (1998). 

M.Salaris and A. Weiss Astron. Astrophys. 335, 943 (1998). 

[3] P. de Bernardis et al. Nature 404, 955 (2000). 

[4] S. Weinberg, Rev. Mod. Phys. 61 (1989) 1. 

[5] A. Guth, Phys. Rev. D 23 (1981) 347; Phys. Lett. 108 B (1982) 389. 

[6] A.D. Linde, Phys. Lett. B 108, 389 (1982); Phys. Lett. B 114, 431 (1982); Phys. 
Lett. B 238 (1990) 160. 

[7] F. Hoyle and J.V. Narlikar, Proc. R. Soc. 273A (1963) 1. 

[8] A.A. Starobinsky, Phys. Lett. B 91 (1980) 99. 

[9] J. Barrow and A.G. Ottewill, J. Phys. A: Math. Gen. 16 (1983) 2757. 

[10] H.-J. Schmidt, Class. Quantum Grav. 7 (1990) 1023; Phys. Rev. D 54 (1996) 7906. 

[11] S. Kluske and H.-J. Schmidt, Astron. Nachr. 317 (1996) 337. 

[12] D. La and P.J. Steinhardt, Phys. Rev. Lett. 62 (1989) 376. 

D. La, P.J. Steinhardt and E.W. Bertschinger, Phys. Lett. B 231 (1989) 231. 

[13] S. Gapozziello, R. de Ritis, G. Rubano, and P. Scudellaro, La Rivista del Nuovo 
Cimento 4 (1996). 



10 



[14] A. Albrecht and P.J. Steinhardt, Phys. Rev. Lett. 48, 1220 (1982). 
[15] R.M.Wald, Phys. Rev. D28 (1983) 2118. 

[16] S.W. Hawking and G.F.R. Ellis, The Large-Scale Structure of Space-Time, Cam- 
bridge Univ. Press (1973) Cambridge. 

[17] S. Capozziello, R. de Ritis, C. Rubano, and P. Scudellaro, Int. Journ. Mod. Phys. D 
4 (1995) 767. 

[18] S. Capozziello, R. de Ritis, and A.A. Marino, Nuovo Cim. 112 B (1997) 1351. 

[19] R.R. Caldwell, R. Dave, P.J. Steinhardt, Phys. Rev. Lett. 80, 1582 (1998). 

[20] G. Magnano, M. Ferraris, and M. Prancaviglia, Gen. Rel. Grav. 19 (1987) 465. 

[21] N.D. Birrell and P.C.W. Davies Quantum Fields in Curved Space (1982) Cambridge 
Univ. Press (Cambridge). 

[22] G. Vilkovisky, Class. Quantum Grav. 9 (1992) 895. 

[23] M. Green, J. Schwarz and E. Witten, Superstring Theory, Cambridge Univ. Press, 
Cambridge (1987). 

A.A. Tseytlin and C. Vafa Nucl. Phys. B 372 (1992) 443. 
G. Veneziano Phys. Lett. B 265 (1991) 287. 

M. Gasperini, J. Maharana and G. Veneziano Phys. Lett. B 272 (1991) 277. 
K.A. Meissner and G. Veneziano Phys. Lett. B 267 (1991) 33. 

[24] P. Teyssandier and P. Tourrenc, J. Math. Phys. 24 (1983) 2793. 

[25] K. Maeda K Phys. Rev. D 37 (1988) 858; Phys. Rev. 39D (1989) 3159. 

[26] D. Wands, Class. Quantum Grav. 11 (1994) 269. 

[27] S. Capozziello, R. de Ritis, A.A. Marino, Class. Quantum Grav. 14 (1997) 3243. 

[28] S. Gottlober, H.-J. Schmidt, and A.A. Starobinsky, Class. Quantum Grav. 7 (1990) 
893. 

[29] A. Vilenkin, Phys. Lett. B117, 25 (1982); Phys. Rev. D30, 509 (1984); Phys. Rev. 
D32, 2511 (1985); Phys. Rev. D33, 3560 (1986); Phys. Rev. D39, 1116 (1989). 

[30] S. Capozziello, R. de Ritis, and A.A. Marino Gen. Rel. Grav. 30 (1998) 1247. 

[31] G. Magnano and L.M. Sokolowski Phys. Rev. D 50 (1994) 5039. 

[32] S. Capozziello and G. Lambiase, Gen. Relativ. Grav. 32 295 (2000). 

11 



[33] I.L. Buchbinder, S.D. Odintsov, I.L. and Shapiro Effective Action in Quantum Grav- 
ity lOP Publishing (1992) Bristol. 

[34] J. Hwang and H. Noh, Phys. Lett. 506 B (2001) 13. 

[35] F.W. Hehl, P. von der Heyde, G.D. Kerhck and J.M. Nester, Rev. Mod. Phys. 48 
(1976) 393. 

[36] A. Trautman Nature 242 (1973) 7. 

[37] V. de Sabbata, Nuovo Cim. A 107 (1994) 363. 

[38] E.W. Kolb, and M.S. Turner The Early Universe 
Addison- Wesley 1990 (Redwood City, Calif.) 

[39] P.J.E. Peebles Principle of Physical Cosmology, 
Princeton Univ. Press 1993 (Princeton). 

[40] V. de Sabbata and C. Sivaram, Astr. and Space Sci. 165 (1990) 51; 
V. de Sabbata and C. Sivaram, Astr. and Space Sci. 176 (1991) 141. 

[41] S. Capozziello and C. Stornaiolo Nuovo Cim. B 113 (1998) 879. 

[42] H. Goenner and F. Miiller-Hoissen, Class. Quantum Grav. 1 (1984) 651. 

[43] S. Capozziello, G. Lambiase, and C. Stornaiolo Ann. Phys. (Leipzig) 10 (2001) 8, 
713. 

[44] M. Tsamparlis Phys. Rev. D 24 (1981) 1451; Phys. Lett. A 75 (1979) 27. 
[45] P. Minkowski, Phys. Lett. B 173 (1986) 247. 



12 



